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ABSTRACT 

We present our measurement of the "bulk flow" using the kinetic Sunyaev-Zel'dovich (kSZ) effect in 
the WMAP 7-year data. As the tracer of peculiar velocities, we use Planck Early Sunyaev-Zel'dovich 
Detected Cluster Catalog and a compilation of X-ray detected galaxy cluster catalogs based on ROSAT 
All-Sky Survey (RASS). We build a full-sky kSZ template and fit it to the WMAP data in W-band. 
Using a Wiener filter we maximize the signal to noise ratio of the kSZ cluster signal in the data. We 
find no significant detection of the bulk flow, and our results are consistent with the ACDM prediction. 
Subject headings: cosmology: cosmic microwave background, cosmology: cosmology: observations - 
(cosmology:) large-scale structure of universe 



I. INTRODUCTION 

Since the flrst claimed de tection of large-scale stream- 
ing by Rubin ct al. ( 1976 ) in Sc galaxies, the issue of 
coherent departures from uniform Hubble flow has been 
the source of much debate. The inflationary model pre- 
dicts that the coherent, large scale peculiar motion of 
matter caused by gravitational potentials, also called 



bulk flow" , is negligible in a ACDM universe ( Strauss & 
Willick|1995 ). This prediction has been tested in the last 
tew decades and it has been the theme of much of the 
work in peculiar velocities. A known flow at small scales ( 
< 60 Mpc h~^) is the motion of the Local Group towards 
a mass concentration of about 10^ Mq known as the 



with the observed CMB dipole (|Kogut et al.|1993 Jarosik 


et al.| 


20111. On large scales. 


Lauer & Postman 


(1994| 



found a non-zero bulk flow at 80/i~^ < R < 110/i^^Mpc 
by using a full-sky peculiar velocity survey consisted of 
119 Abel clusters. But a re-analysis of the data led to 
a reduced bulk flow in a different direction IHudson fcl 



Ebeling ( 1997). And at the same time Riess et al. 
used SNla data and found no evidence of a bulk 
similar scales 



(19971 



how at 



Feldman et al. (2010) used a compilation 
of peculiar velocity surveys and found a non-zero bulk 
flow at i? - lOO/i-^Mpc. 

All methods based on galaxy data are limited by our 
ability to observe and measure accurately at large dis- 
tances. This limits the reach of these methods to scales 
of lOOMpc. An independent measurement method on 
substantially larger scales can shed light on the matter 
and clarify the situation. The largest accessible probe 
to study t he bulk flow to date is the kinetic Sunyaev- 
Zeldovich (Sunyaev & Zcldovich 19801 effect due to the 
coherent motion of clusters ot galaxies with respect to the 
rest frame of the Cosmic Microwave B ackground (CMB) 
radiation (Haehnclt & Tegmark 1996). Peculiar veloci- 
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ties of the electrons in the hot intracluster gas lead to a 
Doppler shift of scattered photons. The shift is propor- 
tional to the product of the line of sight component of 
the peculiar velocity and the electron density integrated 
along the line of sight through the cluster. A coherent 
flow causes an overall dipole seen through various optical 
depths of the clusters. This is a unique pattern that can 
be exploited in measuring the bulk flow. 



Most recently [Kashlinsky et al 
tHe 



( |2008[ ) found a net 

component measured in 



dipole moment m the kSZ 
WMAP data that was consistent wit h a non - zero b ulk 
motion at scales of i? ~ 575/i~^Mpc. 



Kcislcr 



(20091 



peated their analysis and foun d a negligible bulk Eow i n 



contradiction to the findings of Kashlinsky et al. ( 2008 ) 



kyz m WMAP data and again found no significant ve- 



Osborne et al.J^(|2011|) did an independent analysis of the 
locity dipole. 

Although a lot of work has been done on this subject, 
the above summary clearly shows that there is no con- 
sensus on the existence of the bulk fiow, its magnitude, 
depth or direction. In this paper we use two different 
sets of cluster catalogs, based on SZ and X-ray detection 
of the clusters, to measure the bulk flow. Explaining the 
bulk flow within the framework of the ACDM model is 
difficult. Several theoretical models have been proposed 
to explain the possible existence of a large b ulk fiow at 
large scales in the framework of a DGP model ( Wyman & 
Khoury|20l"0l) or by u sing other modified gravi ty models 



( [Afshordi et al.||2ob9[ [Khoury k Wyman 



2009) 



I'he rest of this paper is organized as follows: In Sec- 
tion [2] we briefly describe the data sets used for this anal- 
ysis. In Section [3] we explain the details of the algorithm 
we use in constructing the full-sky kSZ templates, filter- 
ing, template fitting and error estimation. We present 
the results in Section [4] and discuss its implications and 
systematics in Section [5] All through t his paper we 
use th e best-fit ACDM model of WMAP7 ( [Larson et al 
2011) for the cosmological parameters. 



2. DATA 

Our analysis uses two types of data sets: a data map 
and a cluster catalog. We use WMAP data with the 
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Fig. 1. — The spatial distribution of the four catalogs and mask 
used in our analysis. 

highest resolution fuh-sky map to date. And as velocity 
tracers we use two independent catalogs of galaxy clus- 
ters; an SZ selected catalog and an X-ray selected cat- 
alog. Below we explain each of these data sets in more 
detail. 

2.1. CMB Data 

We use co-added inverse-noise weighted data from 
seven single-year maps observed by WMAP at 94 GHz 
(W-band). The maps are foreground cle aned (using the 



foregr ound template model discussed in Hinshaw et al.| 
(l2007[) ) and are at HEALPixFJresolution 9 (TV^idc ^ 512). 
I'he W MAP data arc signa rdominated on large scales, 
/ < 548 ( Larson et al.|201l]) and the detector noise dom- 
inates at smaller scales, ihe noise in WMAP data is a 
non-uniform (anisotropic) white noise that varies from 
pixel to pixel in the map. Pixel noise in each map is 
determined by N^bs with the expression 



(1) 



where ctq is the noise for each differencing assembly and 
can be found on the WMAP data products webpage on 
lambda]^ Nobs is the number of observations at each 
map pixel which is directly proportional to the statistical 
weight, i.e. regions with larger number of observations 
have lower noise variances. Nobs is included in the maps 
available from the LAMBDA website. 

In all of our analysis we use pixel masks to exclude 
foreground-contaminated regions of the sky from the 
analysis. We use a galactic mask which masks 19.30% 
of the sky. We do not mask the point sources. 

2.2. Cluster Catalogs: RASS X-ray Catalogs 

We use two sets of catalogs for the clusters of galaxies: 
three X-ray catalogs and one SZ catalog. There are 627 
clusters in the X-ray catalogs and 189 clusters in the SZ 
catalog. The Clusters in the Zone of Avoidance (CIZA) 
cluster catalog is an X-ray survey that identified clus- 
ters in the region where the magnitude of the galactic 
latitude is less than or equal to 20°. CIZA used X-ray 
data from ROSAT All-Sky Survey (RASS) for its ini- 
tial cluster candidate selection, and accepted or rejected 



|prodi 
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candidate clusters according to optical and near infra- 
red (NIR) observations. This catalog has the locations 
of 73 clusters from the CIZA survey, and the medi an red- 
shift of these clusters is 0.07 ([Ebeling et aL| ( |2002| ). The 
Extended ROAST Brightest Cluster Sample ( eBCS) cat- 
alog is also an X-ray catalog. eBCS was compiled from 
RASS data and identified clusters in the northern hemi- 
sphere with galactic latitudes of magnitude greater than 
or equal to 20°; it is estimated to be 75% complete. Total 
fluxes of clusters in this catalog are between 2.8 x 10~^^ 
and 4.4 x 10~^^ erg cm~^ s~^ (0.1 to 2.4keV). We have 
the locations of 107 clusters from the eBCS cata log, and 
the median redshift of these clusters is 0.13. ( [Ebeling 



et al. (20001). The third X-ray cluster catalog we use is 
tEe~ROMT-ES0 Flux Limited X-Ray (REFLEX) clus- 
ter catalog which covers 4.24 steradians in the southern 
sky. The sample is limited to those clusters with an X-ray 
flux above 3 x 10~^^ erg s~^ cm"^ (0.1 to 2.4 keV), and 
it is estimated to be more than 90% complete. We have 
the locations of 447 clusters identified by RE FLEX, and 
the median redshift of these clusters is 0.09 (Bohringer 
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2.3. Cluster Catalogs: Planck ESZ 

Planck is mapping the whole sky at a few arcminute 
resolution and will eventually produce an SZ detected 
cluster catalog of a few thousand galaxies. What is 
currently available, hereafter Planck ESZ catalog, is an 
isotropic all-sky catalog and has 189 clusters with a me- 
dian redshift of 0.546. Figure [l] shows the spatial dis- 
tribution of the clusters of galaxies in the four catalogs 
discussed above, and Figure [2] compares the stacked clus- 
ter profiles. 

3. GENERATING AND FITTING FULL-SKY KSZ 
TEMPLATES 

We generate a full-sky kSZ template due to the bulk 
flow and measure the components of the bulk flow veloci- 
ties by fitting the template to the data in real space. The 
kSZ template is built using an indicator of the mass of 
the cluster. We don't have a direct mass estimate in the 
catalogs, so we use scaling relations to derive the masses 
based on other properties of the clusters, in particular 
their luminosities and their integrated Compton-K fac- 
tors depending on the catalog. 

3.1. Cluster Parameters: X-Ray Catalogs 

In each of the X-ray catalogs, we have the luminosity 
of each cluster, and we use the luminosity to find the 
virial mass and radius for each cluster. To do so, we first 
find Afsoo = ^[500pc(^)]-R5oqi the mass of the cluster 
contained withm -R500, the radius within which the mean 
overdensity of the cluster is 500 times the critical density 
of the universe, pc{z) = 2.77 5 x W^E'^{z)Mp , Mpc ~^ 
We use the scaling relation of Vikhlinin et al. (2009) to 
compute M500 from the X-ray luminosities 

In Lx = (47.392 ± 0.085) + (1.61 ± 0.14) In Afgoo 
+ (1.850 ± 0.42) \nE{z) - 0.39 ln(/i/0.72) 
±(0.396 ±0.039), (2) 

where Lx is the X-ray luminosity in units erg , M500 
is in units of M0 and E{z) = [17 m(l+z)'^+^Al ^ i s thee x- 
pansion history of the universe ( [Vikhlinin et al. (20091). 
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Fig. 2. — Left) The stacked cluster profiles in Planck ESZ and the X-Ray catalogs in rings of width 0.00145 Radians around each cluster. 
The error bars are given by the standard deviation of this quantity for all the clusters in the catalog. Only clusters with absolute latitude 
\b\ > 25° are used in the stacking. Right) The distribution of M^ir in Planck ESZ and the X-Ray catalogs. 



From M500 we simultaneously calculate the virial mass, 



i.e., mass enclosed within the virial radius, Ry 



and the concentration parameter, c , using the system of 
equations 

47r, 



-[K{z)p,{z)]R 



3 

'vir ' 



(3) 



m^cR^oo/Rvir) 



5.72 



-0.081 



(l + z)0-7i V10"/i-iMq 
where m{x) = ln(l + x) — tt— , c is taken from 

f 



et al 



H-x ' 



Duffy 



2008) based on the N-body simulations with the 
hve-year cosmolo gical parameters and A dz) is 
a function of n„i and 57a ( [Bryan fc Normaii||1998[ ) 



Ac(z) = IStt-" + 82[n{z) - 1] - S9[n{z) 



(4) 



and il{z) — rim(l + z)^/E^{z). Fo r refer ence, 2i?5oo 
approximates Rmr (Komatsu et al. (2011 1). We use 
the concentration parameter to dehne the scaling radius, 
namely 

Ryir /-\ 

rs = —. (5) 

We numerically solve the system of equations ([s]) using 
FuncDesigner library for Pythoi£] 

3.2. Cluster Parameters: SZ Catalog 

Planck ESZ contains the integrated compton-Y at the 
cluster position and within 5i?5oo. Here the compton-Y 
parameter is given by 



P-dL 



(6) 



where (Tt is the Thomson cross section, me is the rest 
mass of an electron, P is the intracluster medium thermal 
electron pressure, and the integral is taken along the line 

^ http:/ /openopt.org/FuncDesignerDoc 



of sight and over the area of the cluster in question. To 
find M500, we use the scaling relation Yji^^g = O.bbYi^fi^gg 
which gives the integrated compton-Y at X-ray position 
and within i?5oo- Furthermore, we use the scaling rela- 
tion 

B 



E{2 



= 10^ 



500 



IQiW, 







(7) 



where A = -4.213, B = 1.72 ([Planck early resuhs. jcl 
(|2011|). The distribution of JVi^,^ in the and X-Kay 



catalogs is shown in Figure [2] Then we use the system 
of equations (|3| to find Myir and Ryir- 

3.3. Cluster Profiles 

We use a /3-model for the cluster profiles. In this model, 
the number density of electrons in a cluster is given by 



ne(r) = rieo 1 + 



2\ -3/9/2 



with 



rieo = 



A-Kr'l (c — arctan c) ' 



(8) 



(9) 



and we take /3 = 2/3 to describe th e X-ray surface bright- 
ness profile of observed clusters (Waizmann & Bartel- 
is defined within the vmal radius ot the 



mann 
cluster 



2009|). 

We calculate Np 



the total number of electrons 
in each cluster, according to 



(10) 



Here Jh = 0.76 is the hydrogen fraction, fgas = 
flf,/Q„i — 0.0168, and rrip is the proton mass. Taking 
the line of sight integrals, the optical depth of a cluster 
in the /3-model is given by 



Tie) 



„ rsneo . _i 
2(Tt — , tan 
' Tip 



\ 



§1 

82 



1 + 



01 

el 



fll) 
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In order to check the effect of the assumed cl uster pro- 
file on our results, we do our analysis using the Navarro, 
Frenk fc White] (|1996[ ) (NFW) model as well. The elec- 
tron density in this model is given by 

n,M = (12) 



1 



with 



n-eo 



47rr3(ln(l + c)-^ 



And the optical depth is given by 



Tie) 



2r.,rteo / 



2rsne 



1 



1 



1 — x- 



tanh 



1+x 



■ tan 



x-l 
l+x 



X < 1 
X = 1 

a; > 1 



where x = ^ (Wright & Brainerd (2000l). The key dif- 
ference between the two models is that the NFW profile 
places most electrons at the central region of the cluster 
whereas the /3-model is more extended throughout the 
cluster and beyond. In both cases we cut off the model 
at Tvir- In practice we project the cluster profiles onto a 
pixelized sphere. To pixelize our models, we calculate 

Ttotai = D^^arNe, (14) 

with ctt the Thomson scattering coefficient. Every clus- 
ter is built within its Ovir- For each pixel in the cluster, 
t{9) is computed and assigned to that pixel. Pixel values 
are normalized in each cluster by 



'^total 



(15) 



Here the sum is taken over all the pixels within the given 
cluster, and flpix is a constant for the entire map given 
47r 

by . Our normalization is such that 



''pix 



'^total 1 



(16) 



where fi is the normalized value of each pixel. The typ- 
ical opitcal depths that we obtain using this algorithm 
are of the order of 10"'^. 

3.4. Template fitting 

The optical depth templates explained above are used to 
create KSZ templates for unit velocities in the a;, y and 
z directions, [1,0,0] km/s, [0,1,0] km/s, [0,0,1] km/s 
respectively by multiplying each pixel by 



cmb 



(17) 



where fi is the vector pointing to the pixel in question, 
and Tcmb is the temperature of the CMB. We also use 
a tem plate for the thermal compo nent of the SZ (tSZ) 
signal 



Sunyaev 



Zeldovich (|1972|) in our analysis 



(18) 



y 



where kB is Boltzman's constant, rrie is the mass of 
an electron, h is Planck's constant and v is the fre- 
quency of the CMB map (94GIIz). We assume a uni- 
form temperature for the electron gas in the clusters to 
be Tg = 4 • l{f K. Before we measure for a bulk velocity, 
we apply a Wiener filter to our templates and data. To 
do so, we take as the denominator of our filter the CMB 
power spectrum plus white noise and we take as numer- 
ator of our filter the power spe ctrum of a kSZ templat e 



(13) with velocity [0,0,1200] km/s ( [Kashlinsky et al.||2008 ) 



Our Wiener filter is based on the theory power spectrum 
and is not derived from the data unlike the Wiener filter 



of Kashlinsky et al. ( 2009). For a s tudy of diffenet filters 



in this context see Osborne et al. (2011) 



Then we use the least square minimization to fit the 



above templates to the data 



X 



(D 



T)*C"i(D- 



T) 



(19) 



where D is the CMB map and a • T = w^rT^*^^^ + VyT^"'' + 

v^T^'^ + AT*^^ is the tSZ plus kSZ fit to the map. We 
take the covariance matrix to be diagonal and determined 
by the pixel noise in WMAP. This is a good approxima- 
tion for the present data because the cluster sample is 
sparse and we are dominated by the pixel noise on small 
angular scales of interest. We use Monte-Carlo simula- 
tions to test our method. The result of the fit is the 
four-vector a that contains the x, y and z components of 
the bulk flow and the tSZ template amplitude. Results 
are presented in Section |4j 

In order to assess the statistical significance of the re- 
sults, we use Monte-Carlo simulations of the CMB sky. 
Simulated maps have three components 

AT(n) - (ATcA/i3(n) + ATtsz{n)) ^ B{n) + N{n), 

(20) 

where ATcmb is a realization of the Gaussian CMB field, 
^rptsz thermal component of the SZ signal simu- 

lated as described in eqn. (18), N{h) is the pixel noise 
and ®B{fi) means convolved with the proper beam of 
the experiment. 

We make 100 realizations of the CMB sky using 
synf ast routine of HEALPisP] with the underlying the- 
ory power spectrum computed with CAMI^ using the 
concordance model. The maps are then convolved with 
WMAP beam for W band. Noise realizations are added 
to the beam convolved maps in the end. Noise maps are 
simulated using eqn. ^ with CTo = 6.549 mK for the W- 
band noise. We test our simulations by comparing their 
average power spectra with the data power spectrum. 

4. RESULTS 

For each cluster catalog we perform the fits using two 
sets of templates: templates based on the NFW model 
for the cluster profiles and templates using the /3-model. 
Using the NFW model, we obtain the bulk velocity [- 
924, 305, -421] ± [864, 852, 644] km/s in the combined 
catalog of X-ray clusters and [-1350, -668, 3520] ± [2581, 
2187, 2085] km/s for the SZ catalog. The results for the 
NFW profile are shown in Fig. |3] Using a /3-model, does 
not significantly change the results. For the /3-model de- 
scribed above we obtain the bulk velocity of [-1110, 123, 

* We use h ealpY, the Python version of HEALPix. 
^ http://cainb.info 
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of the cluster s in the catalog. Ou r results agree with 



a: 



I I X-Ray 

I' j) Planck ESZ 



the results of |Osborne et al.| ([20Tl]); |Kcislcrl (|2009| and 



Velocity Components 

Fig. 3. — Bulk velocity components measured using the two cata- 
logs. Our results do not show a significant bulk flow. The larger er- 
rorbars on the Planck catalog measurements are due to the smaller 
number of objects in the catalog. 

-391] ± [872, 901, 689] km/s for the X-ray catalog and 
[-2940, -1120, 4840] ± [3512, 3443, 2441] km/s for the SZ 
catalog. In order to measure the bulk flow at different 
redshifts, we divided the X-ray catalog into three cata- 
logs at redshifts z < .068, 0.068 < z < 0.13 and z > 0.13. 
The redshift bins are chosen so that all bins have roughly 
the same number of clusters in them. We repeated our 
analysis for each of these samples. The bulk flow was 
consistent with zero in all three redshift bins. 

5. SUMMARY AND CONCLUSION 

We use cluster catalogs and the highest resolution 
maps of the WMAP data (W-band) to measure the bulk 
flow using the kSZ effect. We do not detect a signifi- 
cant bulk flow. The best constraint we get on the bulk 
flow velocities is from the RASS X-ray catalog, assum- 
ing an NFW profile for the clusters. The result is shown 
in Fig. [3| Using a (3 model leads to similar results. 
This is due to the size of the WMAP beam which makes 
the details in the cluster profiles indistinguishable at this 
resolution. The results based on the Planck ESZ cata- 
log have large uncertainties due to the smaller number 



contradict the findings ot Kashlinsky,Atrio-ljarandela fe 



Ebeling (2011 1 where they report a significant detection 
ot a bulk flow inconsistent with zero. It is important to 
have as many independent statistical methods as possible 
to measure and constraint the bulk flow velocities. Fu- 
ture cluster catalogs and CMB maps will help us tighten 
these constraints. 

This analysis can be improved by using larger clus- 
ter catalogs such as the f uture X-ray selected ca talog 
of eRosita All-Sky Survey (Cappelluti et al. 2011), and 
the next releases of Planck SZ selected clusters combined 



with the SZ detected clusters of ACTP ol ( Nieniack et al. 
2010| ) and SPTpol ( [McMahon et al.||2069p . However, 
using more clusters alone will not improve the analysis 
much; one needs to use a higher resolution CMB map 
as well to increase the signal to noise ratio in the mea- 
surement. The only full-sky maps of the CMB at higher 
resolution in the near future are going to be Planck maps. 
Using Planck data with future cluster catalo gs will help 
us pu t tighter constraints on the bulk flow (Mak et al. 
2011 1. As the data get better, it is important to improve 
the theoretical models used for making the templates. 
For example second order effects like the changes in the 
SZ cluster brightness, flux and nu mber counts induced 
by the motion of the Solar System ( Chluba et alj2005t 
and the rela tivistic corrections to the kSZ signal ( INozawa 
et al.|[2006 l need to be taken into account. 
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